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COMPARISON IffiEE EXPERIMENT OF SEVERAL METHODS 
OF PREDICTING THE LIFT OF WINGS IN 
SUBSONIC COMPRESSIBLE FLOW 
By Harry E. Murray 

SUMMARY" 


Several methods of predicting the lift of wings in subsonic cont- 
pressihle flow were compared with experiment . An experimental verifica- 
tion of Eaplan* s formula for the effect of compressibility on the lift 
of wing sections was obtained. 

Semi empirical formulas were developed for predicting the subsonic 
effects of compressibility on the lift of finite— span wings based on 
corrections to the section lift-curve slope. These ssmiempirical 
formulas yielded better agreement with experiment than previously 
derived theoretical methods. The agreement at small sweep angles was 
slightly better when thickness was considered in the semiempirical 
formulas . 

Both experiments and calculations indicated a decrease in the 
variation of lift with Mach number for increasing sweep. 


INTRODUCTION 


The effect of compressibility on the lift of finite— span wings has 
been extensively discussed in previous papers (references 1, 2, and 3, 
for example). These papers discuss compressibility effects in terms 
of an affine transformation based on small-perturbation theory herein 
referred to as the three-dimensional Prandtl transformation. The 
application of the three-dimensional Rrandtl transformation to the 
lifting-line theory of inswept wings is discussed in references 1 and 2. 
An application of Weis singer *s approximate lifting-surface theory of • 
wings of arbitrary sweep (reference k) is discussed in reference 3. 
i/hen compared with experiment, these existing methods did not yield 
entirely satisfactory results. 

Kaplan (reference 5) has shown that inclu ding the thickness of a 
two-dimensional airfoil in calculations of the effect of compressibility 
on the lift results in appreciable effect at high subsonic Mach numbers. 
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Inclusion of thickness in the three-dimensional case consequently may 
also have an appreciable effect. Bather formidable mathematical 
difficulties are encountered in any rigorous attemprb to consider the 
thickness of a finite— span wing' in subsonic compressible flow. To 
attempt an approximate adaptation of Kaplan * s two-dimensional solution 
to finite— span wings therefore seems reasonable. 3n the present paper 
such an approximate adaptation of Kaplan* s results to swept and unswept 
wings is obtained. Tbe method is based on lifting-surface theory (refer- 
ence 6). Available test data are compared with the present method, as 
well as with several other methods which neglect the effect of airfoil 
thickness. 


SKMBOIS 
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finite— span-wing lift coefficient 
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section lift coefficient 
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air mass density 
wing area 
wing chord 
free— stream velocity 
angle of attack 
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section lift-curve slope at M = 0 

section lift-curve Blope at Mach number of M 
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Mach n umb er I 

\Yelocity of sound/ 

sweep angle of wing quarter-chord line, positive for sweepback 
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Ka = 


fa -Me 2 

p = fa - M 2 


^K " 


7 

X 

A 


(1 + |(l - e -2 ^ jx(p. - 1) + |(7 + l)(n 2 - l) 2 J 

ratio of specific teats (1.4 for air) 

airfoil thickness parameter (i}r 0 in reference 7) 

aspect ratio 


A. = A/l - M 2 


E e . lifting— surface-theory correction factor based on A 

E e equivalent value of E Q "based on A c 


A . _i tan A 
A = tan — — — 




[/l - 

sweep factor 


K = 


AHAEZSIS 

The Lift of Two-Dimensional Wings 


Kaplan (reference 5) has presented a method whereby the effect of 
airfoil thickness can he considered in the calculation of the lift of a 
two-dimensional wing at an angle of attack. Kaplan* s formula for the 
effect of compressibility on lift is 


1 





n(n - 1) + ^(7 + D(n 2 - l) 2 J 


( 1 ) 


The effect of thickness is included in this formula by means of the 
parameter X which relates the chord of the airfoil to the radius of 


I 
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the confor mal circle in He odors en* s potential theory of arbitrary 
profiles (reference 7)* As the airfoil thickness approaches zero-, X 
approaches zero and Kaplan’s formula can he seen to approach the Prandtl 
(or Glauert) factor, which is 


1 _ 

F ~ 


ij- = 


1 

1/1 — hi 2 


(2) 


He variations of section lift— curve slope with Mach number for the 
MCA 66,1— 115 airfoil and for the same airfoil with a beveled trailing 
edge (models 1 and 2, table i) are shown, 'in figures 1 and 2, respectively 
(from reference 8). Calculations applying Kaplan* s formula (equation (l)) 
and the Prandtl factor (equation (2)) are also shown for comparison. 

The short vertical lines on the curves of figures 1 and 2 indicate 
the lower limit of the test data below which the experimental curve 
is extrapolated to zero Mach number. Consideration of the thickness 
(the Kaplan method) improves the agreement between theory and experiment . 
Similar agreement is shown in figure 3 for -an airfoil approximating the 
MCA 0012-64 airfoil, designated as E-0009 in reference 9, at zero sweep 
(velocity of free stream normal to wing span) . 


In reference 10, Jcnes indicated that the effect of compressibility 
on a swept two-dimensional wing is the same as an an inswept wing in a 
stream of reduced Mach number M Q , where 

Mq = M cos A (3) 

Kaplan* s formula can thus be adapted directly to the two-dimensional 
swept wing by replacing p by 


^A = 


l/l - Me 2 


l/l — M^cos^a 


(4) 


Calculations made with this modification of Kaplan*s formula, together 
with experimental data and calculations made with the Prandtl 
correction, are presented in figure 3 for the two-dimensional 
MCA 0012—64 airfoil (model 3j table i) at three angles of sweep. 

He results indicate that the foregoing method of predicting the effect 
of compressibility on a two-dimensional swept wing produces good agree- 
ment with experiment for A = 20° when thickness is considered. He 
reason for the somewhat poorer agreement for A = 40° is not known. 

Below the short vertical lines on the curves of figures 1, 2, and 3 
the experimental curve is extrapolated to zero Mach number by means of 
Kaplan *s formula in which p^ is used. 
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Hi© lift of Three-Dimensional Wings 

The effect of compressibility on the lift of three-dimensional -wings 
is- calculated by four different methods. All four methods stem from the 
three-dimensional Prandtl transformation. According to the three-dimensional 
Prandtl transformation as set forth in reference 3, the effect of com- 
pressibility on lift can be obtained by increasing by the factor — — j; 

|/l - M 2 

the lift for. incompressible flow of wings having equivalent aspect ratios 
given by 


Aq = A/l - M 2 


and equivalent sweep angles given by 


tan A c 


tan A 

Vl - M2 


Methods 1 and 2 are based on an interpretation of the three-dimensional 
Prandtl transformation, which is strictly correct only for unswept wings 
of high aspect ratio to which lifting-line theory is applicable. This 
particular interpretation was adapted for application to wings of moderate 
aspect ratio with sweep because it afforded a simple, logical means 
wherein a correction for thickness based on Kaplan* s results (reference 5) 
could be applied. Methods 3 and 4, however, are strict applications of 
the three-dimensional Prandtl transformation. No logical or practical 
method could be discovered for adapting a strict application of the three- 
dimensional Prandtl transformation to the purpose of accounting for the 
effects of wing thickness. 

Method 1 .— According to reference 1, if the three-dimensional 
Prandtl transformation is applied to lifting-line theory, the lift-curve 
slope of unswept wings is 


n - 1 






(5) 


where A c = f3A and |3 = /l — M 2 . Equation ( 5 ) also can be written as 




A + 


fl 57 . 3 ' 
0 * 


(6) 



6 


MCA TH Ho. 1739 


Equation (6) shows, as has already "bean observed in reference 2, that 
for lift ing -line theory the effect of compressibility can he accounted 
for by simply using the correct variation of section lift— curve slope 
with Mach number. Equations (5) and. (6) apply rigorously only to the 
limiting case of unswept wings having very high aspect ratio. Accor ding 
to reference 6, the lift-curve slope of inswept wings in incompressible 
flow can be obtained more correctly than by lifting-line theory from the 
following equation 


Aa± 

' x = AE q + a| 2Li3. 

0 ~ it 


(7) 


which is based, on lifting -lin e theory but corrected according to lifting- 
surface theory. The product AE Q is shown in figure k. If the effect 

of compressibility can be assumed to be accounted for by correcting the 
section lift— curve slope , as in the case of lifting— line theory then 
equation (7) can be written as 


% 



at 57-3 

P * 


AE e +• 


(8) 


Equation (8) is not a strict application of the three-dimensional Prandtl 
transformation because of the presence of the quantity E e . The effect 
of thickness can be approximately introduced into equation (8) as a 
correction to the section lift-curve slope by substituting Pg for ‘P; 
the resulting equation is 





AE 0 


fl 57-3 

Pe 


(9) 


Reference 11 has shown that the effect of sweep can be approximately 
accounted for in formulas for the lift— curve slope of unswept finite- 
span wings by multiplying the section lift— curve slope by the 
factor cos A. If some new factor E A is assumed to account exactly 

for sweep, then equation (9) can be written to Include sweep, as follows: 




57.3 


( 10 ) 


3t 
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For a swept wing the results of the two-dimensional analysis will he 
used to account approximately for the variation of compressibility effect 
with sweep. For this purpose pg- in equation (10) will he based on 

or. In other words, upon the component of Mach number norma], to the 
swept panel. Comparisons of equation (10) with experiment indicated 
that best agreement resulted when A is defined as the sweep of the 
quarter— chord line and 1, as used in the equation for Pg-, corresponds 

to the airfoil section normal to the quarter-chord line. 

Because F^ cannot be accurately evaluated from available theories 
for all useful values of aspect ratio and sweep, equatiarL (10) is 
valuable because of its ability to predict the effect of compressibility 
an a wing for which low-speed experimental data are available. The 
quantity a^F^ can be evaluated from the low^-epeed test condition and 

is assumed to be Independent of Mach number in the subsonic range. 

Equation (10) has been applied as herein explained and is compared 
with experiment in figures 5 to 19 as method 1. The short vertical lines 
on the curves of figures 5 to 19 indicate the lower limit of the test 
data below which the experimental curve is extrapolated to zero Mach 
number by means of method 1. 

Method 2 .— Method 2, which is identical to method 1 except that 
1=0, is included in figures 5 to 19 in order to Indicate the order 
of magnitude of the effect introduced by thickness. When 1=0, Pg- 
is replaced by p in equation (10), and the resulting equation is 




P 


^e + y *a 


2 hi 

it 


Method 3 .— A strict application of the three— dimensional Erandtl 
transformation (reference 3 ) to very thin unswept, or very slightly 
swept, wings transforms equation ( 7 ) to 



Ac^i 


57-3 




( 11 ) 


Equation (11 ) has been applied as a correction to low— speed test data 
for wings of less than ±12° sweep and 1b compared with experiment in 
figures 5 to 11; a^ is evaluated from the low^-speed data by method 1. 
This application is denoted method 3- 
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Method 4 .— The three— dimensional Rrandtl transf ormation can he 
extended to thin, swept wings by means of Welss1nger t s appro xi mate 
lifting-surface theory of reference 4, as explained in reference 3- 
Lift-curve slopes calculated by Weiss1nger*s method can be obtained 
from charts in reference 3 for a wide range of sweep angles, aspect 
ratios, and taper ratios. Because the Weissinger method calculations 
were obtained for a section lift-curve slope of 2n (the thin-alrfoil— 
theory value), Weissinger* s method was not applied as a correction for 
compressibility to low— speed data. The calculated values are compared 
directly with experiment in figures 5 "to 19 and are denoted method k. 


KESUITS AHD DISCUSSIOH 
Comparison of Experiment and Calculation 


Four methods have been used to predict the variation of lift— curve 
slope with Mach number for finite— span wings. Figures 5 to 19 show 
comparisons between these methods and experimental data obtained from 
wind-tunnel tests of wing models. Table I contains supplementary 
information regarding the models and, test conditions for the experimental 
data; this information was obtained from references 8, 9, 12, 13 , 1^, 
and unpublished data. Table I indicates that, for models 4 to 18, jet- 
boundary corrections were either applied or were negligible. All 
corrections were applied in a similar manner to account for the effects 
of boundary— induced velocity and blockage. Failure to apply corrections 
to the data of models 1, 2, and 3 resulted in a Blightly excessive 
increase of lift-curve slope with Mach number at high subsonic Mach 
numbers. Application of the correction would, in general, improve 
slightly the agreement between theory and experiment. 

Examination of figures 5 to 19 shows that all wing tests and all 
four methods of calculation yield an Increase of lift with Mach number. 
Theory and experiment are, therefore, in qualitative agreement. 
Quantitatively, however, the agreement between calculation and experiment 
is not quite consistent. Such Inconsistencies (figs. 8 and 10, for 
example) are to be expected, however, because all the calculations are 
based on potential flow. Methods 1, 2, and 3 ar© applied as corrections 
to low— speed test data which account for the low— speed boundary— layer 
effects. If the Reynolds number changes which usually accompany Mach 
number variation produce variations in the boundary layer of the airfoil, 
an additional variation of lift— curve slope with Mach number, which is 
not predictable by the methods discussed herein, will occur. The marked 
difference between experiment and results calculated by method k in 
figure 6 results from the wing model having a full— span elevator with 
an overhang and an open gap. When the variation of lift-curve slope 
with Mach number was calculated by methods 1, 2, and 3, the gap was 
considered in terms of its effect on the low— speed section lift— curve 
slope. Ho consideration of the gap was made when method 1 was applied. 
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la general, of the four methods, methods 1 and 2 compare most 
favorably with experiment. Of these two methods, method 1 which includes 
a consideration of section thickness appears somewhat better than method 2, 
particularly for small 3 sweep angles. Method 2 gives less increase of 
lift— curve slope with Mach number than method 1. Method 3 indicates 
even less increase of lift— curve slope with Mach number than method 2 
and seems ‘to be in poorer agreement with experiment. Method 4 indicates 
an increase of lift— curve slope with Mach number less than method 2 
but similar to method 3* Coincidence of the calculations of method 4 
and experiment at low Mach numbers was not obtained as was the case with 
methods 1, 2, and 3 because the calculations were al 1 obtained for a 
low— speed section lift— curve slope of 2k. 

The data of figures 5 to 19, either experimental or calculated, 
indicate that the variation of lift— curve slope with Mach number decreases 
as the sweep increases, furthermore, a comparison of methods 1 and 2 
indicates that the effect of thickness upon the variation of lift— curve 
slope with Mach number decreases with sweep, either positive or negative. 


Practical Considerations for Method 1 

• Because the agreement between experiment and method 1 seemed good, 
a practical calculation procedure, based an the assumption that a low- 
speed lift-curve slope is available for the finite— span wing under con- 
sideration, is briefly set forth for this method. The constant a^B^ in 
equation (lb) can be solved for as follows: 


*i*A = = K (12) 

A — dr 

where is the value measured at some low Mach number. The quantity AE e 

is evaluated from figure k; whereas is evaluated from equation (l) 

with pa being used rather than p; that is, only the component of the 
Mach number perpendicular to the wing quarter-chord line is considered. 

The value of K from equation (12) is an effective section lift-curve 
slope for zero Mach number appropriate to the wing 'for which the low- 
speed data were obtained. 

The finite— span lift— curve slope and its variation with Mach number 
for the wing for which K was evaluated can now be written 



AE e 


Eg 

JE 57^3 

* 


(13) 
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If no low— speed lift— curve slope is available from, which to evaluate E, 

Weissinger’s method, which gives C-r at M = 0 (chants of reference 3 ), 

a . 

can "be used. If the' wing is effectively unswept (sweep angles of quarter- 
chord line "between ±12°), E approaches a^, which can "be determ in ed 

from, wing-section tests, many of which are presented in reference 15 . 

The evaluation of fig- requires that the value of X he known for 
the particular airfoil section under consideration. This parameter can, 
of course, he calculated as explained in reference 7 j however, such a 
calculation is often unnecessary because X haB already been calculated, 
as a step in the pressure-distribution calculation, for a large number 
of airfoils. Figure 20 shows the variation of X with thickness for a 
number of airfoil series. The value of X in equation (l) should 
correspond only to the airfoil basic-thickness form (airfoil without 
camber). Equation (l) is, therefore, independent of airfoil camber, 
and figure 20 offers a very wide coverage of airfoils. 

As explained in reference 15 there are two slightly differing groups 
of airfoils both designated as 6— series. The older group shows the low- 
drag range following a comma after the number denoting the chordwise 
position of minimum pressure (model 1, table I, for example). A more 
recent group of airfoils show the low-drag range as a subscript to the 
number denoting the chordwise position of, minimum pressure (model 6, 
table I, for example). The curves of figure 20 apply only to the latter 
group of 6— series airfoils. 


C0NCHJSI0KS 


A conparison with experiment of several methods of predicting the 
lift of wings in subsonic compressible flow indicates the following 
conclusions : 

1. Kaplan* s formula for the effect of compressibility on the lift 
of wing sections is in good agreement with experiment. 

2. Semiempirical formulas derived for finite— span wings agree 
with experiment better than previously derived theoretical -methods. 

3 . A slightly better agreement at small sweep angles results from 
a consideration of thickness in the semiempirical formulas. 
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il-. Both experiment s and calculations indicated a decrease in the 
variation of lift with Mach number for increasing sweep. 


Langley Aeronautical laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Ya., August 6, 19^ 
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figure /.- Vor/af/on of fwo-d/mene/onal I /ft- curve 
c/ope Pv/ff? Mac/? number. Mode/ /. A = 0°. 



d/gure 2.- Var/af/on of iwo - dimensional I iff curve 
s/ope pv/fh Mac/? number . ; Mode/ 2. A=0°. 
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fi/gure 3e Vor/ut/on of two-at/mens/onef / iff- curve 
c/ope with More/? number for three ong/es 
of sweep . Mo c/e/ 3 . 
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F/yure Sr Var/af/oa of f/ff -curve slope with Mach 
numher. A = 0°. A4o7e/ 4. 
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F/p/ure 7- Vernation of tiff curve slope with Mach 
number . A- 0° M octet 6 . 
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F/gc/re /Or Var/af/ 0/7 of //ff -curve dope w/fh Mach 
/7u/7?der. A = 6.3". Mode/ d. 
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figure /S.- Var/af/or? of //f /-curve fiope wrth 
Alacf? rurufer. A. = ~4S.0\ Moc/e/ /4. 
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